Caustic activation of rain showers 
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We show quantitatively how the collision rate of droplets of visible moisture in turbulent air 
increases very abruptly as the intensity of the turbulence passes a threshold, due to the formation 
of fold caustics in their velocity field. The formation of caustics is an activated process, in which a 
measure of the intensity of the turbulence, termed the Stokes number St, is analogous to temperature 
in a chemical reaction: the rate of collision contains a factor exp(— C/St). Our results are relevant 
to the long-standing problem of explaining the rapid onset of rainfall from convecting clouds. Our 
theory does not involve spatial clustering of particles. 
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It is common experience that rainfall can commence 
very abruptly from cumulus clouds (which form when the 
atmosphere is convecting) but has a much slower onset 
from stratiform clouds in a stable atmosphere. This can 
happen even when no part of the cloud is below freezing 
point. It is believed that the difference between convect- 
ing and stable clouds arises because the convection gives 
rise to small-scale turbulent motion which facilitates the 
coalescence of microscopic water droplets ('visible mois- 
ture') into raindrops. This idea has a long history (|l| is 
a significant contribution containing references to other 
early papers), but a satisfying theory has been elusive 
and the topic remains a subject of intensive research, 
reviewed recently in [2j. Numerical experiments have 
shown a dramatic increase in the rate of collision of sus- 
pended particles when the intensity of turbulence exceeds 
a certain threshold. This effect was first described in [3j, 
see also |j] . Here we present a simple quantitative theory 
of this phenomenon, illustrated in figure^ which shows 
the collision rate R of an aerosol as a function of a dimen- 
sionless parameter termed the Stokes number, St, which 
contains information about the turbulence intensity, e, 
and the radius of the water droplets, a (other symbols 
in the caption are defined later). There is a precipitous 
increase in the collision rate at a threshold value of St, 
which was also observed in J, |j] . The current consen- 
sus, represented in |2|, |jj, U, |fj , is that the increased rate 
of collision involves spatial clustering of particles. One 
exception is [(j , which presents a theory having elements 
(described later) in common with our own, but which is 
more complex in its formulation and less precise in its 
conclusions. Spatial clustering plays no role in our the- 
ory. 

The properties of clouds are very variable and the sizes 
of visible moisture droplets have a large dispersion, but 
typically the average of the radius is approximately 10/im 
and the density is n « 10 8 m~ 3 |2j. The motion of the 



droplets is dominated by viscous forces, so that the equa- 
tion for the position r of a droplet is well approximated 
by 

r = j(u(r,t) — r) (1) 

until the particles come into contact 1] (here u(r,t) is 
the velocity field of the air and dots denote derivatives 
with respect to time). According to Stokes's formula for 
the viscous drag on a sphere, 7 = 9p g v/2pfa 2 {v being 
the kinematic viscosity of air, pi and p g the densities of 
water and air respectively): using the values above we 
estimate 7 rs 500 s" 1 . The turbulent motion is a multi- 
scale flow with an approximately power-law spectrum of 
spatial fluctuations of u(r,t) 7\. According to the Kol- 
mogorov theory of turbulence, the short wavelength cut- 
off of the power-law spectrum, 7/ (termed the Kolmogorov 
length) is a function of v and of the rate of dissipation 
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FIG. 1: Particle collision rate as a function of Stokes num- 
ber, for three different values of the Kubo number (defined 
in equation 10). In these simulations we used n = 10 3 , 
a = 2.5 x 1CP 4 , r\ — O.l/V, r = 0.1, and varied Mo and 7: 
Ku = 1 (0), Ku = 0.2 (o), and Ku = 0.04 (□). The theoreti- 
cal curves are (|7J, with actions S obtained from figure^ The 
constant C a in i? a = Cs,na uo/n was fitted (see text). 



per unit mass in the turbulence, e. Similar considerations 
apply to the shortest characteristic timescale, r. Dimen- 
sional arguments [3 then imply 
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The typical velocity difference for points separated by 
the correlation length rj will be denoted by uo- The di- 
mensional arguments of the Kolmogorov theory imply 
that wo ~ ti/t for turbulent flow. The dissipation rate 
e is highly variable between different clouds. Values of 
e = 0.1m 3 s~ 2 are typical |2|; v ~ 10 _5 m 2 s _1 for air in 
the lower atmosphere, giving r\ « 3 x 10~ 4 m, r sa 10 _2 s. 
We can form four independent dimcnsionless param- 
eters describing the microscopic motion of the water 
droplets: these are 
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(d is the dimensionality of space). The values quoted 
above give St sa 0.2. The Kubo number Ku cannot be 
large [2j and the comments above indicate that it is of or- 
der unity for fully developed turbulence, but small values 
of Ku occur in stirred fluids, and this parameter will play 
a role later. From the data above we see that nrf <C 1 
and r\ja >• 1 for water droplets in clouds. 

Saffman and Turner [ij only considered the case of 
St <C 1 in detail, because they believed that St ^> 1 
would only be realised in the most unstable cumulonim- 
bus clouds. Later it was understood that fully-developed 
turbulent motion with high Reynolds number can exhibit 
pronounced intermittency pfl. One consequence is that 
the rate of dissipation can fluctuate wildly relative to the 
mean value, and measurements in clouds have indicated 
that e can have an approximately log-normal distribu- 
tion |2|. The notion that intermittency can promote the 
formation of rain showers by creating localised regions of 
very high turbulent intensity is considered in [2j and [fj. 

Aerosol particles in a turbulent flow may exhibit a de- 
gree of clustering: mechanisms for this process were pro- 
posed in 0, M> la] and reference Q describes recent ad- 
vances. Most of the recent literature on the initiation 
of rainfall by turbulence assumes that the mechanism 
involves clustering [2J. y. |j. |r|. This theory is unsat- 
isfactory for a number of reasons. Firstly, the widely 
accepted view is that clustering is due to particles be- 
ing centrifuged away from regions of high vorticity. It 
is argued that this effect is strongest when St is close 
to unity: it does not occur when the particle motion is 
highly damped (St <C 1), or when the vortices are too 
short-lived (St ^S> 1). Figure ^ shows, however, that the 
collision rate rises abruptly at a threshold value of St 
(which depends upon Ku) and remains high for St ^$> 1: 
this is equally true for simulations of single-scale flows, 
such as figure QJ and multi-scale flows, such as those in 
|3J . Another difficulty is that the clustering can only oc- 
cur for particle separations smaller than 77, but in a cloud 



the density of particles is so low that there is unlikely to 
be more than one particle in a cube of length r). Also, 
it has not been established that the centrifuge mecha- 
nism is sufficiently effective to explain the large increase 
in collision rate. 

In our theory, the dramatic increase in the collision 
rate is a consequence of the formation of 'caustics', illus- 
trated schematically in figure [21 Here the velocity of the 
water droplets as a function of position is initially single- 
valued (curve (a)) but droplets with a large velocity over- 
take slower-moving particles, so that at a later time the 
droplet velocity is multi- valued (curve (b)). The region 
where the velocity is a multi-valued function is bounded 
by two fold caustics. The caustics have two effects which 
could enhance the collision rate of water droplets. Firstly, 
at the caustic lines themselves there can be a divergence 
in the density of particles. This effect is discussed in 
0, llTJl : it is analogous to the divergence of light intensity 
on optical caustics [llj, but we argue that it is not rele- 
vant to collisions between microscopic water droplets in 
clouds, where nrf <§; 1. The other effect of the caustics 
is that when the velocity field is multi- valued, droplets 
at the same position are moving with differing velocity, 
and their relative motion produces collisions. This has 
no analogue in optical caustics. The importance of this 
effect was previously emphasised in [^J. 

Our theory must consider the rate of collision between 
droplets both with and without caustics, and the rate at 
which caustics are formed. Saffman and Turner [ij dis- 
cussed the collision (due to shearing motion) of droplets 
which are advected with the air: this approximation is 
valid in the limit as St — > 0; we term this advective col- 
lision rate R a . At large Stokes numbers a theory due 
to Abrahamson |l2j treats the water droplets as a gas 
of particles with random and uncorrelated motion and 
calculates the collision rate R E by gas-kinetic theory. 
This theory is widely regarded as being applicable only 
at very large Stokes numbers, and indeed the two the- 
ories give results which differ by a factor of order rj/a 
at St s=a 1. This might be taken as evidence that at least 
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FIG. 2: Schematic showing the formation of a caustic: (a) 
Initial configuration, showing manifold representing droplet 
velocity as a function of position, (b) This manifold has devel- 
oped fold caustics. The velocity field is multivalued between 
the folds. 



one of these theories does not work when the Stokes num- 
ber is of order unity. However, we shall argue that the 
gas kinetic model is applicable as soon as caustics have 
formed and that the collision rate is well-approximated 
by R = R^ + f (St, Ku)i? g , where /(St, Ku) is the fraction 
of the coordinate space for which the velocity field has 
become multi- valued due to the formation of caustics. 

First we consider the rate of production of caustics 
in more depth. In earlier works f [lfi Il3j|). two of the 
present authors described the formation of caustics in 
the limit where Ku <C 1. In the one-dimensional case, we 
considered the linearisation of equation QJ, describing 
the small separation in space Sx and velocity Sv for two 
nearby droplets. The equation of motion for X = 5v/5x 
is Q 
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which is a stochastic differential equation, in which the 
velocity gradient acts as a random forcing term. The 
droplets encounter a caustic whenever Sx passes through 
zero, implying that X goes to infinity in one direction, 
then jumps instantaneously to the reflected point at in- 
finity. When fluctuations drive X to a sufficiently neg- 
ative value, it will almost certainly escape to — oo in a 
finite time. In J13j, equation (0J was considered in the 
limit St — > oo, Ku — > 0, where the dynamics of X can 
be approximated by a Langevin process. The associated 
Fokker-Planck equation was solved exactly, giving the 
rate J at which any one droplet passes through caus- 
tics. Similar results were obtained in two [lfj and three 
[9j spatial dimensions. They are expressed in terms of a 
dimensionless parameter X defined in terms of the strain- 
rate correlation function 
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(we use angular brackets to indicate averages). In the 
limit as X — > 0, we find that the rate of caustic formation 
is determined by the rate of escape due to diffusion of the 
coordinate X from an attractive fixed point at X = 0. 
The escape rate is asymptotic to 



J = Jo exp(-5/J) 



(6) 



(for some constant Jo), where in S can be obtained as the 
action of a trajectory of a Hamiltonian function |l5U16j. 
Note that we can write X ~ Ku 2 St, so that JBJ is consis- 
tent with the alternative form quoted in the abstract. In 
one dimension we extract 5=1/6 from our exact expres- 
sion for J, obtained i n Il31 . In two and three dimensions, 
one obtains S sw 0.14 |lfj| and S w 0.12 9] in the limit as 
Ku — > 0. Our numerical results, summarised in figure 
show that the activated escape model for the formation 
of caustics, described by }6"| . is valid even when Ku ap- 
proaches unity. We find that the action S depends upon 



Ku. For our model (described below) we found S = 0.70 
for Ku = 1, S = 0.18 for Ku = 0.2 and S = 0.15 for 
Ku = 0.04, consistent with the limiting value S = 0.14 
for Ku — ► quoted above. 

The dependence of J upon X (or equivalently , on St) is 
analogous to the dependence of the rate of a chemical re- 
action on temperature T, which is well approximated by 
expressions containing an Arrhenius factor, exp(—E/kT) 
(E is the activation energy, k is the Boltzmann constant). 
In this sense we may regard the formation of caustics as 
an activated process, in which the dimensionless inten- 
sity of the turbulence (X or St) plays a role analogous to 
temperature. 

We argue that the function /(St, Ku) describing the 
prevalence of caustics is closely related to the rate of 
caustic formation, so that the collision rate is well ap- 
proximated by 



R = R a + exp(-S/X)R g 



(7) 



Formulae for the advective and gas-kinetic collision rates, 
i? a and R g respectively, are given below. It is equation 
Q which is the theoretical curve illustrated in figure ^ 
using values of S obtained from figure We fitted the 
prefactor of the expression for R a , which is the dominant 
term when St <C 1, but there are no other fitted parame- 
ters. We remark that combining equations (2) and (4) of 
6] gives an expression which can be written in the form 
R = R a + P R c , analogous to our @, in which R c is 
similar to the precise asymptote, R g . The expression for 
the factor P in [6j has a different dependence upon St. 

Our numerical simulations use a two-dimensional 
model in which the divergenceless velocity u(x,t) is 
obtained from a scalar random field ij)(x,t) by writ- 
ing u — (dip/dy,—dip/dx). The field ip(x,t) is a 
Gaussian random function with spatially and tempo- 
rally stationary, isotropic statistics. The mean value is 
zero and the correlation function is (ip(x,t)ip(x' , t')) — 




FIG. 3: Rate of caustic formation as a function of 1 oc Ku 2 St, 
showing evidence for 'activated' behaviour: Ku = 1 (0), with 
action S = 0.70, Ku = 0.2, S = 0.18 (o), and Ku = 0.04, 
S — 0.15 (□). Also shown is the limiting behaviour as Ku 
— > 0, with action S = 0.14 (dashed line). 



Mgexp[ — (t — t') 2 /2t 2 ] exp(— \x — x'\ 2 /2r] 2 ), which gives 
X = ^/7r/2Ku 2 St for this model. The droplets are ini- 
tially randomly positioned, and they are regarded as hav- 
ing collided when their separation falls below 2a (that is, 
we assume that the collision efficiency is unity ij). Fig- 
ure n shows the rate of collision R for a single particle. 
We did not include gravitational settling or other effects 
which also occur without turbulence. The results for our 
single-scale flow are very similar to the simulations using 
Navier-Stokes flows, reported in |i£,|4j, implying that the 
multi-scale aspect of those simulations is not essential to 
their understanding. 

In order to complete the discussion of our formula (0) 
for the collision rate we describe the asymptotic formulae 
for the advective and gas-kinetic collision rates. In the 
advective limit the suspended particles are brought into 
contact by the effect of shearing motion in the flow. The 
typical shear rate is uq/t). The relative velocity for par- 
ticles that will be brought into contact is Av ~ auo/rj, 
and the volume swept in time At is AV ~ a d uoAt/n. 
Collisions typically occur when nAV ~ 1, so that the ex- 
pected rate of collision is R a = C & na d uo/ri (where C a is 
a constant which we fitted numerically). Saffman and 
Turner |l| gave an upper bound on the collision rate 
which can be calculated analytically if correlation func- 
tions of the flow field are known. Their approach gives 
a precise asymptote for the steady hyperbolic flow which 
they discussed, but in general it is an upper bound be- 
cause their calculation does not account for multiple col- 
lisions in generic flows (which are not steady and which 
may be locally elliptic). 

When St ^> I, the inertia of the aerosol particles means 
that their motion becomes uncorrelated with the velocity 
of fluid, and a gas-kinetic model of the type proposed by 
Abrahamson [l2j is appropriate. The maximal Lyapunov 
exponent A for motion of the particles is positive |ljj and 
if the rate of collision satisfies R <C A, then the veloci- 
ties of nearby particles are completely randomised by the 
time they collide (our estimates below verify this assump- 
tion). The aerosol particles become a gas of droplets 
with velocities which are uncorrelated with each other 
and with that of the turbulent air. The velocity distribu- 
tion is Maxwellian, even if u is not Gaussian distributed, 
because the equation of motion for the droplet velocity 
is analogous to an Ornstein-Uhlenbeck process. The rate 
of collision is exactly the same as for a hard-sphere gas 
with the same particle radius and r.m.s. velocity. The 
number of collisions per unit time of a given particle in 
a dilute gas is R g ~ na d ~ 1 v r (R g — 4anv T in two dimen- 
sions), where v T is the mean relative speed of the droplets. 
For a Maxwellian velocity distribution, v T — \J -K j \\j (v 2 ) 
in two dimensions and we have i? g = 2y/irna^/(v 2 ). It 
remains to calculate the r.m.s. velocity of the water 
droplets, (v 2 ), from which we obtain (v 2 ) — 2(v 2 ). The 
solution of the equation of motion (JIJ for the velocity of 
a droplet is 

v(t) =7 / dt' exp[-i(t - t')]u(r{t'),t') . (8) 



The variance of the velocity is obtained by squaring this 
expression and averaging: when St 3> 1 this is asymp- 
totic to 



dt (u(0,t) -u(0,0)) 



(9) 



This gives a collision rate R g ~ na u^y/yf, (for our 
two-dimensional model, i? g = 2 7 ' 4 7r 3 ' 4 nauoi/rr). Thus 
Rg/Rn is of order rj/a when St ss I and Ku « 1. 

In summary, we have shown that the collision rate is 
well approximated by the gas- kinetic model as soon as the 
turbulence intensity X exceeds the action S. Finally, we 
must consider whether this gives a sufficiently rapid rate 
of collision. In three dimensions, the gas-kinetic model 

Qy/^ncPyJlp 2 ). From 10, 



gives a collision rate R g 
we estimate (v 2 ) 



St t]/t. Using the above data, this 
estimate gives a collision rate of i? g ~ 5 x 10 _3 s~ 1 when 
St w 1. We conclude that rainfall can be initiated in a 
timescale of a few minutes, provided that a sufficiently 
large part of the cloud has a turbulence intensity, X which 
exceeds the action for forming caustics, S. 

B. M. acknowledges support by Vetenskapsradet. 



[1 
[2. 
[3. 

K 

[5. 
[6 

[-: 

[8 
[9 

[10 
[11 



[12 
[13. 

[i4; 

[15 
[16] 



P. G. Saffman and J. S. Turner, J. Fluid Mech., 1, 16-30, 

(1956). 

R. A. Shaw, Annu. Rev. Fluid Mech., 35, 183-227, 

(2003). 

S. Sundaram and L. R. Collins, J. Fluid Mech., 335, 75- 

109, (1997). 

Y. Zhou, A. S. Wexler and L.-P. Wang, Phys. Fluids, 10, 

1206-16, (1998). 

M. R. Maxey, J. Fluid Mech., 174, 441, (1987). 

G. Falkovich, A. Fouxon and G. Stepanov, Nature, 419, 

151-154, (2002). 

U. Frisch, Turbulence, Cambridge University Press, 

(1997). 

J. Sommerer and E. Ott, Science, 259, 335, (1993). 

K. Duncan, B. Mehlig, S. Ostlund and M. Wilkinson, 

Phys. Rev. Lett, 95, 240602, (2005). 

M. Wilkinson and B. Mehlig, Europhysics Lett, 71, 186- 

92, (2005). 

M. V. Berry, Singularities in Waves, Les Houches Lecture 

Series Session XXXV, eds. R. Balian, M. Kleman and J- 

P Poirier, North Holland: Amsterdam, 453-543, (1981). 

J. Abrahamson, Chem. Eng. Set., 30, 1371-9, (1975). 

M. Wilkinson and B. Mehlig, Phys. Rev. E, 68, 040101, 

(2003).bbbbb 

L. Piterbarg, SIAM J. Appl. Math., 62, 777-800, (2001). 

B. Mehlig and M. Wilkinson, Phys. Rev. Lett., 92, 

250602, (2004). 

M. I. Freidlin and A. D. Wentzell, Random Perturbations 

of Dynamical Systems, Springer, New York, (1984). 



